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We study a certain class of classical one dimensional piecewise linear maps. For these systems we 
introduce an infinite family of Markov partitions into equal cells. The symbolic dynamics generated 
by these systems is described by bistochastic (doubly stochastic) matrices. We analyze the structure 
of graphs generated from the corresponding symbolic dynamics. We demonstrate that the spectra 
' of quantized graphs corresponding to the regular classical systems have locally Poissonian statis- 

, tics, while quantized graphs derived from classically chaotic systems display statistical properties 

characteristic of Circular Unitary Ensemble, even though the corresponding unitary matrices are 
sparse. 
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■ I. INTRODUCTION 
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Although graphs served as models of physical systems for a long time, the interest in properties of spectra of the 
^\ • corresponding quantum systems has rather short history. Quantum graphs were introduced as a "toy" model for studies 
Y \ of quantum chaos by Kottos and Smilansky They observed that spectral statistics of fully connected graphs is 

well reproduced by random matrix theory (RMT) and may be explained using the exact trace formula. Schanz and 
Smilansky |^] proved that the spectral correlation function of simple graphs coincides with CUE expression for 2x2 
|Jh matrices. The trace formula served to calculate the spectral correlation function and to explain the deviations from 
Ph the RMT predictions for star graphs [Q. Tanner introduced directed graphs and the corresponding unitary transfer 
matrices and demonstrated fast convergence of their spectral properties to RMT statistics with increasing matrix 

■ size. A model for study of the Anderson localization by means of quantum graphs was introduced by Schanz and 
J> ' Smilansky while other properties of quantum graphs were recently analyzed in 0,^. Continuous dynamics on 

graphs was recently discussed by Barra and Gaspard g| , who generated Markov process by means of Poincare surface 
1^ of section. 

•—^ . In this work we propose another approach leading to quantum graphs. We are interested in a certain class of one 
I ' dimensional (ID) piecewise linear maps with infinite family of Markov partition into equal cells. For these systems we 
, define a family of the corresponding graphs by means of bistochastic (doubly stochastic) transfer matrices B. If the 
transfer matrices are unistochastic (i.e. there exists a unitary matrix U such that Bij — \Uij\'^) the graphs may be 
\ quantized. We may thus establish a family of the unitary transfer matrices by varying the lengths of bonds of graphs 
. a-nd investigate their spectra. We demonstrate a link between the character of dynamics of the classical system and 

the spectral properties of the corresponding ensemble of unitary matrices. For a class of regular ID systems we show 
that the statistics of the spectrum is locally Poissonian. If the classical system has an arbitrary small, positive KS- 
^ ^ entropy, the spectral statistics of the corresponding ensemble of unitary matrices displays fluctuations characteristic 
^> , to Circular Unitary Ensemble (CUE) |l^Jl^. These matrices are sparse (only few nonzero elements in each row and 
^ ' column) and by construction differ from typical random matrices generated according to the Haar measure on U (N). 
C3 ; The paper is organized as follows. In section || we introduce the class of ID maps of our interest. Section |l| 
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contains an analysis of the structure and spectral properties of quantized graphs corresponding to regular systems. In 
section we study graphs corresponding to Bernoulli shift and to some systems with arbitrary small, but positive 
metric entropy. Section ^ is devoted to study of generic chaotic systems, for which we find CUE like spectral 



properties of the corresponding ensembles of unitary matrices. The concluding remarks are presented in section VI 
while appendix |^ contains the condition necessary for a bistochastic matrix to be unistochastic and the calculation of 
free parameters for averaging over different graphs with the same topology. In appendix |b] we show the construction 
of unitary matrices corresponding to the classical map studied in section |V|. 

II. ONE DIMENSIONAL DYNAMICAL SYSTEMS 

Consider a discrete one dimensional mapping acting on an interval I d M 

Xn+l = f{Xn), (1) 
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where / : / ^ / is piecewise linear. Assume that / fulfils three conditions: 

(i) there exists a Markov partition of the interval / on M equal cells Ei,i — 1 . . . M and / is linear on each cell Ei , 

(ii) for any y (z I 

1 



E 



^e/-i({y}) 



= 1, 



(2) 



(iii) the finite transfer matrix B, describing the evolution of measures uniform on the cells of Markov partition, under 
the action of the system / (the Frobenius-Perron operator) is unistochastic. 

Condition (i) means that the graph of the function / (later called diagram of / to avoid the confusion with graphs 
corresponding to the map /) is composed of straight lines (neither horizontal nor vertical) with endpoints onto M x M 
grid on the square I x I. Some examples of functions fulfilling this condition are plotted in figure |l]. If we rescale / 



to [0, 1] and write f{x) — CjX + bj for x € Ij where IJ 



I, then all cj have to be nonzero integers and also, the 
numbers bj and must be rational (multiplication by a sufficiently large M makes all these numbers integer). 

Any subpartition of the interval / into N — Mk equal cells (fc S IN) is also a Markov partition, so we have an infinite 
family of Markov partitions for our models. The Frobenius-Perron operator, reduced to densities (measures) constant 
on each cell Ej, may be described by a matrix B of finite size M. Let us write the elements of B explicitly. The 
system fulfilling the condition (i) is fully characterized by a set of M linear functions 



fi{x) = CiX + bi 
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(3) 



By construction and Mbi are both integers. It is convenient to express and bi in terms of two integers li and 
as 



li. bi = 



TJii — i{mi — li) 
M ' 



It is now easily seen that the only non-zero elements of B in the i-th row are 

1 



Bi 



\mi -k\' 



min(m.i, U) + 1 < j < max(m.i, li). 



(4) 



(5) 



Closer inspection of the structure of the above formula reveals the connection between the structure (positions of the 
non-zero elements) of the matrix B and the diagram of the system / rotated clockwise by 7r/2. 

The element Bij is equal to the probability of going from the cell Ei to Ej. From (H) B is stochastic (X^^i ^ij — l)i 
all probabilities of going out of the cell Ei sum to unity. The probability of coming to the cell Ej from Ei is equal to 
mean l/|/'(a;)| over x e Ei and f{x) G Ej. Thus due to condition (ii) the probabilities of transitions to any fixed cell 
Ej also sum to one (X)il=i ^ij — l)i the transition matrix B is bistochastic. It is well known that all eigenvalues of 
a bistochastic matrix are located in the unit circle, and the leading eigenvalue is equal to one |12|. The condition (iii) 
will be used in section 




FIG. 1. Examples of piecewise linear functions defining classical maps with Markov partition consisting of (a) M = 6, (b) 
M = 2 and (c) M = 6 equal cells. Only system (a) fulfils (i)-(iii). Observe that the ID system shown in panel (b) is regular 
even though the average slope exceeds unity, this system does not fulfil the condition (ii), the matrix B corresponding to the 
system (c) is not unistochastic, so (iii) is violated. 
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We may compute the KS-entropy of Markov chain generated by transition matrix B |13| ] 

N N 

Hks{B)^-Y, J2 ' (6) 

1=1 j=i 

where p is the normahzed left eigenvector of B corresponding to the leading eigenvalue, pB = p, ^^^iPi = 1- In 
condition (i) we require that / is linear (not only piecewise linear) on each cell Ej. Then the Markov partition on 
M equal cells is a generating partition of the system, so Eq. (^) gives the dynamical entropy of the system. The 
transition matrix is bistochastic in our case, so all components of p are nonzero {pi — 1/N) and Hks = if and only 
if all Bij e {0, 1}. It is easy to see that if Vx G / : = 1 the system is regular. We have showed that this is also 

a necessary condition, provided the conditions (i) and (ii) are fulfilled, since only the systems with |/'| = 1 have all 
Bij g {0, 1}. Note that in general there exist regular systems with slopes larger than one, but they do not satisfy the 
condition (ii). An example of such a system is plotted in Fig. |l|(b). 



III. REGULAR SYSTEMS 



Partition of the interval / determines the transition matrix Bij describing the symbolic dynamics of the ID map f{x). 
We introduce a directed graph with vertices corresponding to each cell Ei of the partition and bonds corresponding 
to all nonzero elements of the transition matrix Bij . Changing the number N of the cells in the partition we change 
the dimension of the matrix B^-^\ so we may obtain a family of graphs corresponding to any ID map fulfilling the 
condition (i). 
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FIG. 2. An example of a regular ID system defined by function / plotted in left top diagram, the corresponding graphs and 
the transition matrices representing Frobenius-Perron operators for Markov partitions into A = 3 and A = 9 cells. The dots 
represent zero elements of matrices. Note a link between the position of the nonzero elements of the matrices S'^' and the 
shape of the diagram of / (rotated clockwise by 7r/2). 



If we consider the regular systems with |/'| ~ 1 the only nonzero elements of matrices Bij are equal to 1. The 
matrices of this type contain only one nonzero element equal to 1 in each row and column and act as a permutation 
on a vector with A^ elements. The corresponding graph has only one incoming and one outgoing bond at each vertex, 
so its valency is equal to one. It is composed from the closed loops. We say that there is an inversion at vertex j 
of the graph if /' = — 1 on the corresponding cell Ej of the partition. Dividing the interval / into M equal cells we 
find that the longest loop of the graph has at most M vertices. If we consider a subpartition of / into A^ = Mk cells, 
each vertex of the graph will be replaced by fc vertices of a new graph. The loops with an even number of inversions 
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will form k new disconnected loops of the same period. The loops with an odd number of inversions will create [k/2] 
disconnected loops with period doubled ([•] denotes the integer part of a number) and also, if k is odd, one loop of 
the initial period. Fig. |2| shows an example of a regular system, for which the minimal Markov partition contains 
M = 3 cells: Ei, E2, E3. Due to /' = — 1 in Ei there exists one inversion in the 3-vertices graph corresponding to 
this system. For larger N the graph is composed of [k/2] loops with 6- vertices and one 3-vertices loop for an odd k. 
The classical map has one periodic orbit with period 3 and a continuous infinite family of orbits with period 6. Note 
similarity in the spectrum of periods of periodic orbits of the map and the graph for large odd k. 

The quantization of a directed graph may be related to problem of finding a unitary matrix U, such as |f/ij P = Bij 
for all indices i,j The matrix U may be seen as a Floquet operator governing the discrete time evolution, while t/" 
represents the evolution operator at discrete time t = n. Analyzing quantum maps, corresponding to time dependent 
periodical systems, one studies the statistical properties of phases of the unimodular eigenvalues of the unitary matrix 
U The elements of U may be expressed by Uij — ry e*^'^ , where the modulus are fixed by the transition 

matrix — y/Bij, and the phases Lij , corresponding to lengths of the bonds of the graph, are not specified by the 
corresponding ID map. We only impose that they fulfil the condition of unitarity U'^U — 1. Varying the phases Lij 
we obtain a family of graphs - an ensemble of unitary matrices. This corresponds to varying the lengths of all bonds 
of the graph with its topology preserved. We are interested in spectral properties of the ensemble of unitary matrices 
defined in this way. If S is a permutation matrix, there is only one nonzero element in each row and each column of 
U, so all phases Lij appearing in matrix U are uncorrelated. 

The N xN unitary matrix U corresponding to a regular ID system is composed of M square diagonal or antidiagonal 
blocks of size k. For the system presented in Fig. || the matrix U has the form 

/ Vi\ 

U= \ V2 \ , (7) 
\0 V3 J 

where M = 3, (Vi)™; = 5m,k+i-ie''^"' is an antidiagonal k x k matrix, (^2)™/ = Smie'-^'"+'° and (Va)™/ = Smie''^"'+^'' 
are diagonal matrices. There are N free phases Lm (the column index is superfluous here so we omit it). The 
eigenvalues of block matrix U may be expressed by means of the spectra of its blocks Vi [|l5| 



eigj^nv^^j, (8) 

where the M-th root is a M valued function. Unitary matrix V — YifLi ^ of size k is diagonal or antidiagonal, 
depending on parity of the number of inversions in the corresponding loop in the basic M-vertex graph. The nonzero 
elements of V are equal to e*^^ , where Cj, j — I . . . k are sums of different phases Lm, so Cj are uncorrelated. If V is 
diagonal the spectrum of U consists of ^e^{^i+'^™)/^'iy-^- The eigenphases of U belonging to the same interval 
of the length 27r/M stem from different series Cj, so they are uncorrelated and the spectrum has locally Poissonian 

statistics. If V is antidiagonal, its spectrum contains V e^^^i^^>'+'^-^^ and the spectrum of U is composed of roots of 
unimodular complex numbers of order 2M . The character of the spectrum remains Poissonian on the intervals of the 
length tt/M. 



eig(t/) = ^^yd^ 



M 

\ 




FIG. 3. Spectral statistics for regular systems: (a) level spacing distribution P(s), (b) spectral ridigity A3(L) and (c) spectral 
form factor K{t). The average is taken over an ensemble of 10^ matrices of the size A*' = 150 with phases (bond lengths) chosen 
randomly. The statistics reveal deviations from the predictions of the Poissonian ensemble (dot dashed line) due to correlations 
on larger length scales. 

The above argument may be generalized for the situation, in which the permutation matrix B is reducible, composed 
of several permutations with different periods Mi . . . Mp. In this case the eigenphases of U are uncorrelated in intervals 
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of length Tr/maxMi. Fig. ^ shows the commonly used spectral statistics calculated for the system shown in 

Fig. y. Note that the spectral ridigity A3(L), which measures correlations in the spectrum at the length scale L, 
saturates at L w 25 = iV/6 (there are N = 150 eigenphases rescaled to the mean density constant and equal to 1). 
The divisor 6 is equal twice the number of cells in minimal Markov partition, since an inversion in one of the cells 
occurs {f'{x) = —1 iov X € El). 



IV. SIMPLE CHAOTIC SYSTEMS 



The Bernoulli shift, /(x) = 2x mod 1, is one of the simplest ID maps displaying chaotic dynamics. This system 
fulfils the conditions (i), (ii) and (iii). The graphs corresponding to this system were analyzed by Tanner The 
system, graphs and matrices are plotted in Fig. Existence of a certain interval for which |/'(x)| > causes the mean 
valency of the graphs to become larger than 1. Some vertices of the graphs have more than one outgoing bond, and 
this fact is related with a positive topological entropy of the corresponding classical map. The number of such vertices 
grows proportionally with N and their presence is reflected in statistical properties of spectra of the corresponding 
unitary matrices. For the binary graphs all vertices have two outgoing bonds and as shown by Tanner |^] the spectral 
statistics converges fast to the predictions of CUE. The traces of the unitary evolution operator [/"■ may be expressed 
by a sum over periodic orbits of the corresponding classical graph, analogous to the semiclassical Gutzwiller trace 
formula p^ ]. For quantum graphs the formula 

Tr L/" = A^e^" , (9) 

v£VO(n) 

is an identity, where VOin) is the set of periodic orbits of period n, the amplitude = Y\(ij)i=u \f^i ^'^^ length 
is the sum of lengths of bonds belonging to the periodic orbit v. 
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FIG. 4. Bernoulli shift, the corresponding binary graphs and the matrices representing Frobenius- Perron operators for Markov 
partitions into N = 2 and A*' = 4 cells. Dots represent matrix elements equal to zero. For larger the structure of B^^' 
reflects the shape of f{x). 



We will show that every periodic orbit of the graph corresponds one to one to a periodic orbit of the dynamical 
system. Assume that v = jij2 • • • is a periodic orbit of period n of the graph, where ji denote the vertices belonging 
to v and |/'| > 1 on at least one of the cells Ej^. Let the functions 

Fjr- E,,^, 3 X ^ f-\x) e E,^ (10) 

be the local inverse of /, where we identify ji and The function Fj^ is a contracting mapping since < i. 

We define 

F = Fj, oF,,o...o (11) 
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which is a contraction of the cell Ej-^ onto itself, because | < 1 for % ^ c. The function F corresponds to moving 
along the orbit backwards. Due to the Banach theorem there exists exactly one fixed point of F in E^^ , namely 
Xjj = F{x^). Hence there is a periodic orbit of period n in the dynamical system corresponding to the orbit v on the 
graph. The inverse is also true, so there is one to one correspondence between the orbits of the ID chaotic system 
and the corresponding graphs. 

The sum ^ goes over periodic orbits of the classical chaotic map. The amplitudes A^, built from probabilities Bij, 
are equal to Hr^To^ so they are equal to the inverse square roots of the instabilities of the orbits of 

the system. Only the lengths are not determined yet, so we can average over the lengths of the bonds imposing the 
constraints due to unitarity of quantum propagator U . With growing dimension the spectral statistics of quantized 
binary graphs tends rapidly to the CUE statistics. Tanner observed it analyzing spectral form factors for graphs with 
rationally independent bond lengths Fig. |^ shows the spectral statistics averaged over an ensemble of lO'^ matrices 
of iV = 100, corresponding to graphs with randomly chosen bond lengths (see appendix ^). 




FIG. 5. Spectral statistics for Tanner binary graphs averaged over bond lengths: (a) level spacing distribution P(s), (b) 
spectral ridigity A3(L) and (c) spectral form factor K{t). The CUE prediction are represented by dot dashed lines (which 
almost coincide with the numerical data). 
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FIG. 6. Chaotic system with Hks = 5 ln2 and M = 2, the function / defining the system, the corresponding graphs and 



the transition matrices plotted for Markov partition into A*' : 
dots. 



4 and N = 8 cells. Zero elements of the matrices are marked by 



We may also analyze a system, for which |/'(a;)| > 1 only in certain regions Ic C /. The KS-entropy of such a 
system is equal to the entropy of a system equivalent to f{x) in Ic times the ratio vol(/c)/vol(/). This allows us 
to construct a system with an arbitrary small, positive KS-entropy. Consider the system presented in Fig. ^ with 
dynamical entropy Hks = 5 1'^2 (equal to the topological entropy), as an illustration of such situation. The structure 
of the unitary evolution operator U is the same as the structure of the bistochastic transition matrix B which is 
determined by the shape of the diagram of /. Block structure visible in Fig. allows us to compute the spectra as 
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in Eq. (|8|), where V is the product of all blocks of U. If there is an even number of antidiagonal blocks, the resulting 
matrix V is a product of one matrix with Poisson distribution of the spectrum and a certain number of non-diagonal 
matrices with a CUE-like spectrum. In ref. we analyzed the spectra of composed ensembles consisting of products 
of random matrices, each pertaining to a given ensemble of unitary matrices. The CUE-like spectrum was found to 
be robust with respect to the multiplication by matrices of other ensembles. This is related to the fact that the Haar 
measure (CUE) is invariant with respect to multiplication. Following these lines of arguments, we conclude that the 
spectrum of V has a CUE-like properties. The spectrum of U is equal to the Af-th root of the spectrum of V (M 
denotes the number of blocks in U) 



e 



i't'u ^^^(i',+2i^)/M ^ j = l.../c, Z = 0...M-1, (12) 



where e^"^^ is the j-th eigenvalue oiV , it consists of M copies of the spectrum of V contracted M times and placed 
one after another, so its statistical properties are locally like these of CUE. Fig.J^ shows the statistical properties of an 
exemplary chaotic system presented in Fig. ^, for which all spectra consist of M = 2 identical parts of length tt. Note 
that the level spacing distribution and the spectral ridigity (for small L) show good agreement with the predictions of 
RMT. On the other hand the spectral form factor (i.e. the Fourier transform of the two point correlation function of 
the spectral density ||l^) exhibits huge fluctuations since it measures correlations in the entire spectrum and reveals 
the presence of 2 copies of the same spectrum. However, the average taken over a small interval oi t = n/N displays 
a CUE-like behaviour. 




FIG. 7. Spectral statistics for quantum graphs corresponding to the system shown in Fig. |6|: (a) level spacing distribution 
P(s), (b) spectral rigidity A3(L) and (c) spectral form factor K{t) averaged over 10^ matrices of size A'' — 200. Dot dashed 
lines represent CUE statistics. The form factor displays large fluctuations, although the average over a window of r of size 
At = 0.07 (solid line) is close to CUE results. 



Local agreement with the predictions of RMT is preserved even though there exist an inversion in the map /, which 
implies an antidiagonal blocks in B^^\ To show this we analyzed a system consisting of 4 blocks, presented in Fig. ||, 
while the corresponding results are shown in Fig. ||. 
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FIG. 8. Chaotic system with Hks = j ln2 and M — A, with an antidiagonal blocks in S''^': function / defining the system, 
corresponding graph, and transition matrix. Find a link between its structure and the shape of f{x). 



In the argument presented in this section we have assumed that both parts of the system (one regular with \f' {x) \ = 1 
and the other chaotic with |/'(a;)| > 1) are connected, so the probability of going from one to the other is positive. If 
this is not true, the system splits into two subsystems, one regular and one chaotic. The spectrum of such a system 
will contain a mixture of eigenvalues of both subsystems, so the level spacing statistics could be described by an 
appropriately adopted Berry- Robnik distribution [ p^ . 
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FIG. 9. As in Fig. Q for the system shown in Fig. ^. 



V. GENERIC SYSTEMS 



A generic ID piecewise linear system fulfilling the conditions (i), (ii) and (iii) is chaotic and |/'(x)| > 2 for all x 
in /. The mean valency of the corresponding graphs is also equal to or greater than 2 and the topological entropy of 
these graphs -fftop ^^t smaller than In 2. The topological entropy of the map / is equal to the log of the spectral 

radius (absolute value of the largest eigenvalue) of the connectivity matrix B [^|, where Bij = 1 if B^j > 0, and 
Bij = otherwise. It is worth to consider, which of these graphs one may quantize by constructing corresponding 
family of unitary matrices. It is equivalent to ask, which bistochastic transition matrices i3, are unistochastic, so 
there exist unitary U, such that \Uij\'^ = Bij. It is clear that not all dynamical systems fulfilling (i) and (ii) give rise 
to bistochastic Frobenius-Perron matrices, which are unistochastic. As an example let us consider the system with 
M = 6 shown in Fig. |](c) 

fi{x) = 6x, for < a; < i, 

/2,3(x) =3a;- i, fori<a;<|, (13) 
f'i,5A^) = 2x - 1, for 2 < 2; < 1. 

The corresponding transfer matrix 

111111" 

2 2 2 
2 2 2 

3 3 ^ ' 
3 3 
3 3 . 

is clearly bistochastic but not unistochastic. Indeed one checks that the necessary conditions for unistochasticity are 
not fulfilled (see Appendix ^ Eq. (|A(j), fc = 3 and / = 5). This is the smallest matrix corresponding to a dynamical 
system of the considered type with this property, some other may be constructed for M = 28 = 14 + 7 + 4 + 2+1 or 
other perfect numbers. In appendix |^ we present a necessary condition for unistochasticity and discuss some features 
of the ensemble of matrices U corresponding to the same bistochastic matrix B. 
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FIG. 10. Classically chaotic ID map ("four legs map"), corresponding graph and the transition matrices for Markov partitions 
into = 4 and N = S cells. The slope of the map / in each cell Ei determines the size of elements in each row of the matrix 
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One of the maps studied by Dellnitz et al. [|19| and showed in Fig. |lO|, may serve as an example of a classical 
map, for which transition matrices B are unistochastic (see Appendix ^). Its KS-entropy is equal to | In 2, while the 

topological entropy equals to In 3. Instead of computing the spectral radius of B'-^^ it is enough to observe that each 
point in / has exactly three preimages, and the topological entropy may be computed as log of the mean number of 
preimages pO| . The average should be taken with respect to the Parry measure of maximal entropy jl^ , which is not 
easy to specify in the general case. However, since the number of preimages is constant, its knowledge is not necessary 
in this case to obtain the required result i/top = In 3. Numerical investigations performed for the corresponding 
ensemble of 10^ matrices of size N — 200 confirmed a CUE-like character of spectral statistics in the entire circle of 
eigenvalues as shown in Fig. pi]. 




FIG. 11. As in Fig. W for the four legs map shown in Fig. hj 



Instead of averaging over an entire ensemble of unitary matrices with random phases, we may study spectral 
statistics for one concrete member of the ensemble. Figure 12 shows the spectral statistics for a single matrix U, in 
which the phases, not determined by unitarity, being set randomly (see Appendix ^). Observe an agreement with 
CUE predictions (the form factor K{t) is averaged over the window of r. At = 0.08). Interestingly, for this system a 
CUE-like behaviour was also found for a single matrix Uq (corresponding to the map plotted in Fig. |l^ with all free 
parameters set to zero (thus exhibiting strong correlations between phases). Matrix Uq is real (therefore orthogonal). 
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FIG. 12. As in Fig. ^ for the four legs map shown in Fig. |lo[ without averaging over the ensemble of unitary matrices. The 
statistics calculated from single matrix of size 2500 with randomly chosen free phases are plotted. 



Although a single member of the ensemble displays CUE-like spectrum, the averaging over free phases provides 
larger statistics. The form factor may be expressed as a double sum over periodic orbits 

if(r) = l(|TrC/"|2) = l/ J2 ^.^.'e^^^-^^-'A, (15) 

\ e-po(n) / 

where r = n/N and the averaging goes over all allowed realizations of the graph (all bond lengths Lij, for which 
unitarity of U is preserved). The difference of the lengths of orbits is equal to — L^r = J2(ij)eu ~ 'l2{im)£u' ^im- 
Inserting new variables L,j = — Ln — Lij + Ln we have 

Ly — Li,' — {Lij + Lii + Lij — Lii) — ^ (i;„j + L;i + Li„j — ill), (16) 

where the lengths Lij, i,j = 2...N are fixed by the condition assuring unitarity of U (see Appendix while 
{2N — 1) remaining phases Ln and Lij are drawn randomly with a uniform distribution on [0, 27r). Due to averaging 
over Lii, Lij the only non vanishing terms in the sum (^5|) arise from periodic orbits v and z/' of the same degeneracy 
class which visit each vertex of the graph the same number of times. This observation generalizes results of Tanner 
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for binary graphs Therefore the calculation of K{t) can be reduced to a combinatorial problem after finding 
the fixed lengths . The calculations of form factor for fully connected graphs j^] and binary graphs j|] shows fast 
convergence of the spectral statistics to RMT predictions with increasing matrix size N. 

VI. CONCLUSIONS 

We proposed to study a certain class of piecewise linear ID dynamical systems, for which the Markov transition 
matrix (of size M) is bistochastic. It has a structure of the diagram of /, defining the classical map, rotated clockwise 
by 7r/2. For each subpartition of the Markov partition we construct a graph and a family of unitary matrices, of size 
equal to a multiplicity of M . This procedure might be thus considered as a "quantization" of a ID classical maps, in 
a sense that it produces a family of corresponding unitary matrices. The above quotation marks represent the fact 
that no canonical action-angle quantization of classical systems with ID phase space is possible. 

We demonstrated that the spectral properties of ensembles of unitary matrices defined in this way reflect the 
character of the classical dynamics. For regular maps the graphs are composed of several finite disconnected loops 
and the spectral statistics of the corresponding unitary evolution operators is (locally) Poissonian. On the other 
hand, chaotic systems (with arbitrarily small, but positive metric entropy), lead to ensembles of unitary matrices with 
(locally) CUE- like spectral fluctuations. These unitary matrices are sparse (only a few non zero elements in each row 
and column), in contrast to the random matrices, typical with respect to the Haar measure on U{N). 

We have shown that the classical ID system and the corresponding family of graphs are topologically conjugated, 
since each periodic orbit of the map corresponds to exactly one periodic orbit of each graph. Applying finer and finer 
Markov partitions to the system, introduces a natural limit in the family of the corresponding graphs. It is the limit 
of number of vertices tending to infinity with the set of periodic orbits preserved. This seems to be the correct limit to 
reach random matrix theory predictions for the spectral statistics of quantized graphs. Spectral properties of a graph 
may be expressed by the exact trace formula based on a sum over periodic orbits. This sum goes over all periodic 
orbits of the ID system and the amplitudes are equal to inverse square roots of the instabilities of each orbit. The 
phases in the sum remain free parameters of the model, with constraints due to unitarity of the evolution operator. We 
find conditions necessary for a bistochastic transfer matrix to be unistochastic, which allow us to perform averaging 
over random phases (bond lengths) with unitarity preserved. 
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NOTE ADDED IN PROOF 

After this work was completed we learnt about very recent related works on unistochastic matrices of Berkolaiko |2^] 
and Tanner H]. 

APPENDIX A: UNISTOCHASTIC MATRICES 

For a real bistochastic matrix B of size N 

N N 

Bjk>0, Y^Bjk^l, J2B,k^l, (Al) 
j=i fe=i 

we want to find a unitary matrix U, such that iC/jfcp = Bjk- We can manipulate with N'^ variable phases Ljk of the 
matrix U 

Ujk = y/B~^ e'"^^" . (A2) 

The unitarity of U requires that 
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N 

E 



Jkl 



(A3) 



Equations (A3) for k — I are fulfilled due to stochasticity of B, Eg . (Rl^ ). Equations (A3) for k < I are complex 
conjugated of these for k > I. Transforming the system of equations (|A3D 



uiu^'u - i)t/^ = uu\uu'^ - 1) = 



C/C/^ -1 = 



(A4) 



we found that {N — 1) equations of the system are depending on the other ones, because the diag onal equations of 
the righthand side are equivalent to conditions for bistochasticity of i? - the first sum in Eq. (|Al| ). We subtracted 
the unity from N since there is one common condition stemming from each of the sums in Eq. ( |Al| ), X]iJ=i ^ij ~ ^ ■ 
The number of independent real equations in (A2) equals thus {N — 1)^. 

We introduce new variables L^j = Lkj — Lki — Lij + Ln, so the equations (A3) are equivalent to 



N 

E 



(A5) 



Since L^i = iij = 0, so the {N — 1) equations ( |A3| ) depend on (A^ — 1) real independent variables Lkj, k, j — 2 . . . N . 
These equations give the values of L^j ■ For N = 2 there exists only one phase L22 equal to tt as shown by Tanner . 
This value does not depend on the bistochastic matrix Bjk- Eqs. ( [A5| ) do not depend on L^i and Lij, so we may 
average over them, assuming the uniform distribution over the unit circle, P{L) = ^ for L € [0, 2tt). This is equivalent 
to averaging of the matrix U ~ D1UD2 over diagonal unitary random matrices Di and D-?, with a fixed matrix U, 
fulfilling \U,k\''^B,k 



The system of equations ( A3 ) is solvable only if 



max V nh 

n=l...N 



V BkmBl, 



1 



N 



< 2 2^ V^kjBi^ and 



max v B„ 

m=l...N 



i" Br, 



1 ^ 
- 2 E V^j fc 



(A6) 



which is the condition for the possibilit y of balance the maximal term in the sums (A3) by the rest and the analog 
condition for the sums stemming from (A4). This necessary condition for a bistochastic matrix to be unistochastic 
becomes sufficient in 3 x 3 case. For = 2 all bistochastic matrices are orthostochastic (and therefore also unis- 
tochastic). However for A^ > 3 sev eral examples of bistochastic matrices, which are not unistochastic, are known 
they do not satisfy the condition ([A^). Statistical properties of random unistochastic matrices are studied in 



APPENDIX B: CONSTRUCTION OF A FAMILY OF QUANTUM MAPS 

The construction of quantum graphs presented in section ^ hinge on the possibility of finding an unitary (or 
orthogonal) matrix corresponding to the bistochastic matrix obtained for a given dynamical system for an arbitrary 
Markov subdivision of the original partition. Rather then attempting the formulation of the most general theorem 
concerning the existence of such unistochastic matrices, we present here a construction for a particular case of the 
"four legs map" presented in Fig. |lO|. In other words we show that the all bistochastic matrices B^^'^ corresponding 
to this classical dynamical system are unistochastic. 

The original bistochastic matrix for the "four legs map" introduced in p!9| for the initial Markov partition into four 
cells reads 



2 2 

B(4) -1(1111 
4 1 1 1 1 1 
2 2 



(Bl) 



Let us consider the following choice of a corresponding unitary matrix 



^(4) 





fV2 


-V2 





' \ 


1 


1 


1 


1 


1 


2 


1 


1 


-1 


-1 




\ 





V2 


-V2 J 



(B2) 
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The bistochastic matrix for the two times finer subdivision reads 



(22 

1 



4 



\ 
2 2 
1 



1 1 

1 1 1 1 

1 1 1 1 



1 
2 
Vo 



1 1 

2 



1 


2 2/ 



(B3) 



Using our knowledge of the unitarity of we want to find a unitary matrix U^^\ such that b\-' = I^Yj This 
can be done by rephcating twice the rows of C/(4) with an appropriate shift. Thus U'^^^ may read as follows 



t/(8) 



I V2 


^V2 

















\ 








^/2 


-V2 














1 


1 


1 


1 


























1 


1 


1 


1 


1 


1 


-1 


-1 


























1 


1 


-1 


-1 














V2 


-x/2 








V 

















V2 


-V2/ 



(B4) 



Note that the number of the non zero elements in C/*^®^ is two times larger than the number of the non zero elements 
in U^'^\ By replicating each row of the matrix U^'^^ k times we obtain in the similar way unitary matrices of the size 
4fc. 

The success of this construction is due to the fact that the nonzero elements of the first and fourth rows of C/W are 
equal, which implies orthogonality of the following pairs of rows in (2,3), (2,5), (4,7) and (6,7). However this 

is not the only way of constructing the family of unitary matrices U^^\ For example we found another family t/'-^' 
originating from 



[7(4) = t 

2 



( V2 -\/2 
1 
1 








1 -1 
1 1 -1 

V2 V2 ) 



(B5) 



even though the procedure described above fails in this case. This fact allows us to hope that the class of classical ID 
maps leading to quantum graphs (fulfilling the conditions (i), (ii) and (iii)) is not trivial. 
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